Cavity quantum electrodynamics for superconducting electrical circuits: an 

architecture for quantum computation 
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We propose a realizable architecture using one-dimensional transmission line resonators to reach 
the strong coupling limit of cavity quantum electrodynamics in superconducting electrical circuits. 
The vacuum Rabi frequency for the coupling of cavity photons to quantized excitations of an adjacent 
electrical circuit (qubit) can easily exceed the damping rates of both the cavity and the qubit. 
This architecture is attractive both as a macroscopic analog of atomic physics experiments and 
for quantum computing and control, since it provides strong inhibition of spontaneous emission, 
potentially leading to greatly enhanced qubit lifetimes, allows high-fidelity quantum non-demolition 
measurements of the state of multiple qubits, and has a natural mechanism for entanglement of 
qubits separated by centimeter distances. In addition it would allow production of microwave 
photon states of fundamental importance for quantum communication. 

PACS numbers: 03.67.Lx, 73.23. Hk, 74.50.+r, 32.80.-t 



I. INTRODUCTION 



Cavity quantum electrodynamics (cQED) studies the 
properties of atoms coupled to discrete photon modes 
in high Q cavities. Such systems are of great interest 
in the study of the fundamental quantum mechanics of 
open systems, the engineering of quantum states and 
the study of measurement-induced decoherence 0, 0, Q , 
and have also been proposed as possible candidates for 
use in quantum information processing and transmission 
0, H Q. Ideas for novel cQED analogs using nano- 
mechanical resonators have recently been suggested by 
Schwab and collaborators 0,0 • We present here a realis- 
tic proposal for cQED via Cooper pair boxes coupled to a 
one-dimensional (ID) transmission line resonator, within 
a simple circuit that can be fabricated on a single micro- 
electronic chip. As we discuss, ID cavities offer a number 
of practical advantages in reaching the strong coupling 
limit of cQED over previous proposals usingdiscrete LC 
circuits @ B, large Josephson junctions 0, 0, E3, or 
3D cavities (lTL liH IT^ . Besides the potential for en- 
tangling qubits to realize two-qubit gates addressed in 
those works, in the present work we show that the cQED 
approach also gives strong and controllable isolation of 
the qubits from the electromagnetic environment, per- 
mits high fidelity quantum non-demolition (QND) read- 
out of multiple qubits, and can produce states of mi- 
crowave photon fields suitable for quantum communica- 
tion. The proposed circuits therefore provide a simple 
and efficient architecture for solid-state quantum com- 
putation, in addition to opening up a new avenue for 
the study of entanglement and quantum measurement 
physics with macroscopic objects. We will frame our dis- 
cussion in a way that makes contact between the lan- 
guage of atomic physics and that of electrical engineer- 
ing. 



We begin in Sec. [H] with a brief general overview of 



cQED before turning to a discussion of our proposed 
solid-state realization of cavity QED in Sec. Hill We then 
discuss in Sec.|W]the case where the cavity and the qubit 
are tuned in resonance and in Sec.[V]the case of large de- 
tuning which leads to lifetime enhancement of the qubit. 
In Sec. IVII a quantum non-demolition read-out protocol 
is presented. Realization of one-qubit logical operations 
is discussed in Sec. IVIII and two-qubit entanglement in 
Sec. I Villi We show in Sec. |D[] how to take advantage 
of encoded universality and decoherence-free subspace in 
this system. 



II. BRIEF REVIEW OF CAVITY QED 

Cavity QED studies the interaction between atoms 
and the quantized electromagnetic modes inside a cav- 
ity. In the optical version of cQED 2], schematically 
shown in Fig. ^a), one drives the cavity with a laser 
and monitors changes in the cavity transmission resulting 
from coupling to atoms falling through the cavity. One 
can also monitor the spontaneous emission of the atoms 
into transverse modes not confined by the cavity. It is 
not generally possible to directly determine the state of 
the atoms after they have passed through the cavity be- 
cause the spontaneous emission lifetime is on the scale of 
nanoseconds. One can, however, infer information about 
the state of the atoms inside the cavity from real-time 
monitoring of the cavity optical transmission. 

In the microwave version of cQED Q , one uses a very 
high Q superconducting 3D resonator to couple photons 
to transitions in Rydberg atoms. Here one does not di- 
rectly monitor the state of the photons, but is able to 
determine with high efficiency the state of the atoms af- 
ter they have passed through the cavity (since the excited 
state lifetime is of order 30 ms). From this state-selective 
detection one can infer information about the state of the 
photons in the cavity. 
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FIG. 1: (color online), a) Standard representation of cavity 
quantum electrodynamic system, comprising a single mode of 
the electromagnetic field in a cavity with decay rate k, cou- 
pled with a coupling strength g — £ rrns d/h to a two-level 
system with spontaneous decay rate 7 and cavity transit time 
(transit, b) Energy spectrum of the uncoupled (left and right) 
and dressed (center) atom-photon states in the case of zero 
detuning. The degeneracy of the two-dimensional manifolds 
of states with n — 1 quanta is lifted by 2g\/n + 1. c) Energy 
spectrum in the dispersive regime (long dash lines). To sec- 
ond order in g, the level separation is independent of n, but 
depends on the state of the atom. 



The key parameters describing a cQED system (see 
TableHJ are the cavity resonance frequency w r , the atomic 
transition frequency f2, and the strength of the atom- 
photon coupling g appearing in the Jaynes-Cummings 
Hamiltonian nil 



H = Hu> T 
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Here H K describes the coupling of the cavity to the con- 
tinuum which produces the cavity decay rate k = oj t /Q, 
while H 1 describes the coupling of the atom to modes 
other than the cavity mode which cause the excited state 
to decay at rate 7 (and possibly also produce additional 
dephasing effects). An additional important parameter 
in the atomic case is the transit time (transit of the atom 
through the cavity. 

In the absence of damping, exact diagonalization of the 
Jaynes-Cumming Hamiltonian yields the excited eigen- 
states (dressed states) [l5j 



-,n) 



cos0 n |J.,n) + sin0 n |T,n+l) (2) 
-sin0 n |J.,n) + cos0 n |T,n + l) (3) 
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and A = Q — w r the atom-cavity detuning. 

Figure ^Jd) shows the spectrum of these dressed-states 
for the case of zero detuning, A = 0, between the atom 
and the cavity. In this situation, degeneracy of the pair 
of states with n quanta is lifted by 1g\Jn + 1 due to the 
atom-photon interaction. In the manifold with a sin- 
gle excitation, Eqs. J5J) and © reduce to the maximally 
entangled atom-field states |±~0) = (||, 1) ± ||,0)) /y/2. 
An initial zero-photon excited atom state ||, 0) will there- 
fore flop into a photon ||, 1) and back again at the vac- 
uum Rabi frequency g/ir. Since the excitation is half 
atom and half photon, the decay rate of |±, 0) is (k+7)/2. 
The pair of states |±, 0) will be resolved in a transmission 
experiment if the splitting 2g is larger than this linewidth. 
The value of g = £ rms d/h is determined by the transition 
dipole moment d and the rms zero-point electric field 
of the cavity mode. Strong coupling is achieved when 
g > k, 7 [Hj. 

For large detuning, g/A 1, expansion of Eq. (@J 
yields the dispersive spectrum shown in Fig.^i). In this 
situation, the eigenstates of the one excitation manifold 
take the form 15] 

|—) ~ -( 5 /A)U,0) + |T,l) (7) 
|+^) - II, 0) + (g/A) |t, 1). (8) 

The corresponding decays rates are then simply given by 



T^y ~ (, 9 /A) 2 7 + At 

~ 7 + ( 5 /A) 2 At. 



(9) 
(10) 



More insight into the dispersive regime is gained by 
making the unitary transformation 



U = exp —(aa + — a'cr ) 



(11) 



and expanding to second order in g (neglecting damping 
for the moment) to obtain 



UHU^ w h 
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As is clear from this expression, the atom transition is ac- 
Stark/Lamb shifted by (g 2 /A)(n + 1/2). Alternatively, 
one can interpret the ac-Stark shift as a dispersive shift 
of the cavity transition by a z g 2 / A. In other words, the 
atom pulls the cavity frequency by ±g 2 / kA. 
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parameter 


symbol 


3D optical 


3D microwave 


ID circuit 


resonance/transition frequency 


lo t /2tt, Q/2ty 


350 THz 


51 GHz 


10 GHz 


vacuum Rabi frequency 


g/n, g/ujr 


220 MHz, 3 x 10~ 7 


47 kHz, 1 x 10 -7 


100 MHz, 5 x 10~ 3 


transition dipole 


d/eao 


~ 1 


1 x 10 3 


2 x 10 4 


cavity lifetime 


1/k,Q 


10 ns, 3 x 10 7 


1ms, 3 x 10 8 


160 ns, 10 4 


atom lifetime 


l/ 7 


61 ns 


30 ms 


2 /ns 


atom transit time 


^transit 


> 50 fis 


100 fjs 


CO 


critical atom number 


No = 2jK/g 2 


6 x 10" 3 


3 x 10" 6 


< 6 x 10" 5 


critical photon number 


mo = 7 2 /2ff 2 


3 x 10~ 4 


3 x 10" 8 


< 1 x 10~ 6 


# of vacuum Rabi flops 


"Rabi = 2g/ (k + 7) 


~ 10 


~ 5 


~ 10 2 



TABLE I: Key rates and cQED parameters for optical Q and microwave Q atomic systems using 3D cavities, compared 
against the proposed approach using superconducting circuits, showing the possibility for attaining the strong cavity QED 
limit (riRabi S> 1)- For the ID superconducting system, a full-wave (L = A) resonator, uj T /2n — 10 GHz, a relatively low Q of 
10 4 and coupling j3 = C g /C E = 0.1 are assumed. For the 3D microwave case, the number of Rabi flops is limited by the transit 
time. For the ID circuit case, the intrinsic Cooper-pair box decay rate is unknown; a conservative value equal to the current 
experimental upper bound 7 < 1/(2 [is) is assumed. 



III. CIRCUIT IMPLEMENTATION OF CAVITY 

QED 

We now consider the proposed realization of cavity 
QED using superconducing circuits shown in Fig. [21 A 
ID transmission line resonator consisting of a full-wave 
section of superconducting coplanar waveguide plays the 
role of the cavity and a superconducting qubit plays the 
role of the atom. A number of superconducting quantum 
circuits could function as artificial atom, but for definite- 
ness we focus here on the Cooper pair box [|| 0, 0, 0] . 



A. Cavity: coplanar stripline resonator 

An important advantage of this approach is that the 
zero-point energy is distributed over a very small effec- 
tive volume (w 10~ 5 cubic wavelengths) for our choice 
of a quasi-one-dimensional transmission line 'cavity.' As 
shown in appendix ^ this leads to significant rms volt- 
ages V® ms ~ \Jhjj-cjcL between the center conductor and 
the adjacent ground plane at the antinodal positions, 
where L is the resonator length and c is the capaci- 
tance per unit length of the transmission line. At a res- 
onant frequency of 10 GHz {hv/ks ~ 0.5 K) and for a 
10 fim gap between the center conductor and the adja- 
cent ground plane, V rms ~ 2 /iV corresponding to elec- 
tric fields £ rms ~ 0.2 V/m, some 100 times larger than 
achieved in the 3D cavity described in Ref. [3|- Thus, 
this geometry might also be useful for coupling to Ryd- 
berg atoms [l9| . 

In addition to the small effective volume, and the fact 
that the on-chip realization of cQED shown in Fig. [21 
can be fabricated with existing lithographic techniques, a 
transmission-line resonator geometry offers other practi- 
cal advantages over lumped LC circuits or current-biased 
large Josephson junctions. The qubit can be placed 
within the cavity formed by the transmission line to 




FIG. 2: (color online). Schematic layout and equivalent 
lumped circuit representation of proposed implementation of 
cavity QED using superconducting circuits. The ID transmis- 
sion line resonator consists of a full-wave section of supercon- 
ducting coplanar waveguide, which may be lithographically 
fabricated using conventional optical lithography. A Cooper- 
pair box qubit is placed between the superconducting lines, 
and is capacitively coupled to the center trace at a maxi- 
mum of the voltage standing wave, yielding a strong electric 
dipole interaction between the qubit and a single photon in 
the cavity. The box consists of two small (~ 100 nm x 100 nm) 
Josephson junctions, configured in a ~ 1 /nm loop to permit 
tuning of the effective Josephson energy by an external flux 
$ext- Input and output signals are coupled to the resonator, 
via the capacitive gaps in the center line, from 50 Q trans- 
mission lines which allow measurements of the amplitude and 
phase of the cavity transmission, and the introduction of dc 
and rf pulses to manipulate the qubit states. Multiple qubits 
(not shown) can be similarly placed at different antinodes 
of the standing wave to generate entanglement and two-bit 
quantum gates across distances of several millimeters. 



strongly suppress the spontaneous emission, in contrast 
to a lumped LC circuit, where without additional special 
filtering, radiation and parasitic resonances may be in- 
duced in the wiring |2(ij| . Since the resonant frequency of 
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the transmission line is determined primarily by a fixed 
geometry, its reproducibility and immunity to 1 /f noise 
should be superior to Josephson junction plasma oscilla- 
tors. Finally, transmission line resonances in coplanar 
waveguides with Q ~ 10 6 have already been demon- 
strated [2ll . l22l | , suggesting that the internal losses can 
be very low. The optimal choice of the resonator Q in 
this approach is strongly dependent on the intrinsic decay 
rates of superconducting qubits which as described be- 
low, are presently unknown, but can be determined with 
the setup proposed here. Here we assume the conserva- 
tive case of an overcoupled resonator with a Q ~ 10 4 , 
which is preferable for the first experiments. 



B. Artificial atom: the Cooper pair box 

Our choice of 'atom', the Cooper pair box is a 

mesoscopic superconducting island. As shown in Fig. |3 
the island is connected to a large reservoir through a 
Josephson junction with Josephson energy Ej and ca- 
pacitance Cj. It is voltage biased from a lead having 
capacitance C g to the island. If the superconducting gap 
is larger than both the charging energy E c = e 2 /2Cs 
(where Cs = Cj + C g is the total box capacitance) and 
temperature, the only relevant degree of freedom is the 
number of Cooper pairs N on the island. In this basis, 
the Hamiltonian describing the superconducting island 
takes the form 

H Q = 4E c J2(N-N g ) 2 \N)(N\ 



Hamiltonian then reduces to a 2 x 2 matrix 



N 



Ej 



^(|A+1) (N\+h.c.) 



(13) 



N 



where N g — C g V g /2e is the dimensionless gate charge 
representing the total polarization charge injected into 
the island by the voltage source. 

In the charge regime, 4E C ^> Ej, and restricting the 
gate charge to the range N g 6 [0, 1], only a pair of ad- 
jacent charge states on the island are relevant and the 



Ej.Cj 



FIG. 3: Circuit diagram of the Cooper pair box. The gate 
voltage is connected to the island through an environmental 
impedance Z(ui). 



H, 



E„ 



Ej 



Q 



(14) 



with E e i — 4i?c(l — 2N g ). The Cooper pair box can in 
this case be mapped to a pseudo spin-1/2 particle, with 
effective fields in the x and z directions. 

Replacing the Josephson junction by a pair of junctions 
in parallel each with energy Ej/2, the effective field in 
the x direction becomes Ej cos(7r$ cxt /<f> )/2. By thread- 
ing a flux <& ext in the loop formed by the pair of junctions 
and changing the gate voltage V g , it is possible to con- 
trol the effective field acting on the qubit. In the setup 
of Fig. [3 application of dc gate voltage on the island can 
be conveniently achieved by applying a bias voltage to 
the center conductor of the transmission line. The res- 
onator coupling capacitance Co , the gate capacitance C g 
(the capacitance between the center conductor of the res- 
onator and the island) and the capacitance to ground of 
the resonator then act as a voltage divider. 



C. Combined system: superconducting cavity QED 

For a superconducting island fabricated inside a res- 
onator, in addition to a dc part V^ dc , the gate voltage 
has a quantum part v. As shown in appendix 1X1 if the 
qubit is placed in the center of the resonator, this latter 
contribution is given by v = V® ms (a' + a). Taking into 
account both V$ c and v in 114(1. we obtain 



He 



-2E C (1 



2nf)a* 



2N g -a z ). (15) 



Working in the eigenbasis I IT) , ||)} of the first two terms 
of the above expression |23( , and adding the Hamiltonian 
of the oscillator mode coupled to the qubit, the Hamilto- 
nian of the interacting qubit and resonator system takes 
the form 



H = haj T ( a* a + 



1 




(16) 



2N g - cos(9)a z + sm(9)a x ). 



Here, a x and a z are Pauli matrices in the eigenbasis 
{|t), ||)}, 9 = &Tctan[Ej/4E c (l ~ 2A g dc )] is the mix- 
ing angle and the energy splitting of the qubit is £1 = 

Note that contrary to the 
2 



[4£ c (l-2iV g dc )] 2 |2i 

case of a qubit fabricated outside the cavity where the Ng 
term in 1(13(1 has no effect, here this term slightly renor- 
malize the cavity frequency w r and displaces the oscillator 
coordinate. These effects are implicit in Eq. ((16(1 . 

At the charge degeneracy point (where N g = 
C e Vg C /2e = 1/2 and 9 = tt/2), neglecting rapidly os- 
cillating terms and omitting damping for the moment, 
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Eq. (|16fl reduces to the Jaynes-Cummings Hamiltonian 
l[TJl with Q = Ej and the vacuum Rabi frequency 

9 = T\I^l> (17) 

where [3 = C g /C^. The quantum electrical circuit of 
Fig. |2 is therefore mapped to the problem of a two-level 
atom inside a cavity. Away from the degeneracy point, 
this mapping can still be performed, but with a coupling 
strength reduced by sin 9 and an additional term propor- 
tional to (a* + a). 

In this circuit, the 'atom' is highly polarizable at the 
charge degeneracy point, having transition dipole mo- 
ment d = hg/£ lms ~ 2 x 10 4 atomic units (eao), or 
more than an order of magnitude larger than even a 
typical Rydberg atom 15]. An experimentally realis- 
tic [l8| coupling j3 ~ 0.1 leads to a vacuum Rabi rate 
g/ir ~ 100 MHz, which is three orders of magnitude 
larger than in corresponding atomic microwave cQED 
experiments 0, or approximately 1% of the transition 
frequency. Unlike the usual cQED case, these artificial 
'atoms' remain at fixed positions indefinitely and so do 
not suffer from the problem that the coupling g varies 
with position in the cavity. 

A comparison of the experimental parameters for im- 
plementations of cavity QED with optical and microwave 
atomic systems and for the proposed implementation 
with superconducting circuits is presented in Table [I] 
We assume here a relatively low Q — 10 4 and a worst 
case estimate, consistent with the bound set by previous 
experiments with superconducting qubits (discussed fur- 
ther below), for the intrinsic qubit lifetime of I/7 > 2 /is. 

The standard figures of merit |24j for strong coupling 
are the critical photon number needed to saturate the 
atom on resonance tuq = j 2 /2g 2 < 1 x 10 -6 and the 
minimum atom number detectable by measurement of 
the cavity output A^ = 2^n/g 2 < 6 x 10~ 5 . These re- 
markably low values are clearly very favorable, and show 
that superconducting circuits could access the interesting 
regime of very strong coupling. 



IV. ZERO DETUNING 
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FIG. 4: Expected transmission spectrum of the resonator in 
the absence (broken line) and presence (full line) of a super- 
conducting qubit biased at its degeneracy point. Parameters 
are those presented in Table The splitting exceeds the line 
width by two orders of magnitude. 



in the proposed superconducting implementation is so 
strong that, even for the low Q = 10 4 we have assumed, 
2<7/(k + 7) ~ 100 vacuum Rabi oscillations are possi- 
ble. Moreover, as shown in Fig. the frequency split- 
ting (g/n ~ 100 MHz) will be readily resolvable in the 
transmission spectrum of the resonator. This spectrum, 
calculated here following Ref . [2{| , can be observed in the 
same manner as employed in optical atomic experiments, 
with a continuous wave measurement at low drive, and 
will be of practical use to find the dc gate voltage needed 
to tune the box into resonance with the cavity. 

Of more fundamental importance than this simple 
avoided level crossing however, is the fact that the Rabi 
splitting scales with the square root of the photon num- 
ber, making the level spacing anharmonic. This should 
cause a number of novel non-linear effects [T^j to appear 
in the spectrum at higher drive powers when the average 
photon number in the cavity is large ((n) > 1). 

A conservative estimate of the noise energy for a 10 
GHz cryogenic high electron mobility (HEMT) ampli- 
fier is riamp = kBTfqlhio ~ 100 photons, where Tjv is 
the noise temperature of the amplification circuit. As a 
result, these spectral features should be readily observ- 
able in a measurement time t meas = 2n amp /(n)n, or only 
~ 32 /is for (n) ~ 1. 



In the case of a low Q cavity (g < k) and zero de- 
tuning, the radiative decay rate of the qubit into the 
transmission line becomes strongly enhanced by a factor 
of Q relative to the rate in the absence of the cavity . 
This is due to the resonant enhancement of the density 
of states at the atomic transition frequency. In electrical 
engineering language, the ~ 50 fl external transmission 
line impedance is transformed on resonance to a high 
value which is better matched to extract energy from the 
qubit. 

For strong coupling g > k, 7, the first excited state be- 
comes a doublet with line width (k + j)/2, as explained 
in section [H] As can be seen from Table U the coupling 



V. LARGE DETUNING: LIFETIME 
ENHANCEMENT 

For qubits not inside a cavity, fluctuation of the gate 
voltage acting on the qubit is an important source of 
relaxation and dephasing. As shown in Fig. [21 in prac- 
tice the qubit 's gate is connected to the voltage source 
through external wiring having, at the typical microwave 
transition frequency of the qubit, a real impedance of 
value close to the impedance of free space (~ 50 f2). 
The relaxation rate expected from purely quantum fluc- 
tuations across this impedance (spontaneous emission) 
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(18) 



where SV(+f2) = 2M7 Re[Z(fi)] is the spectral density of 
voltage fluctuations across the environmental impedance 
(in the quantum limit). It is difficult in most experiments 
to precisely determine the real part of the high frequency 
environmental impedance presented by the leads con- 
nected to the qubit, but reasonable estimates 0] yield 
values of T\ in the range of 1 /is. 

For qubits fabricated inside a cavity, the noise across 
the environmental impedance does not couple directly to 
the qubit, but only indirectly through the cavity. For 
the case of strong detuning, coupling of the qubit to the 
continuum is therefore substantially reduced. One can 
view the effect of the detuned resonator as filtering out 
the vacuum noise at the qubit transition frequency or, in 
electrical engineering terms, as providing an impedance 
transformation which strongly reduces the real part of 
the environmental impedance seen by the qubit. 

Solving for the normal modes of the resonator and 
transmission lines, including an input impedance R at 
each end of the resonator, the spectrum of voltage fluc- 
tuations as seen by the qubit fabricated in the center of 
the resonator can be shown to be well approximated by 



S V {Q) 



2hu> T k/2 



Lc A 2 



{K/2f 



(19) 



Using this transformed spectral density in (|18fl and as- 
suming a large detuning between the cavity and the 
qubit, the relaxation rate due to vacuum fluctuations 
takes a form that reduces to 1/Ti = j K = (g/A) 2 K ~ 
1/(64 /is), at the qubit's degeneracy point. This is the re- 
sult already obtained in Eq. (fTTTfl using the dressed state 
picture for the coupled atom and cavity, except for the 
additional factor 7 reflecting loss of energy to modes out- 
side of the cavity. For large detuning, damping due to 
spontaneous emission can be much less than k. 

One of the important motivations for this cQED ex- 
periment is to determine the various contributions to the 
qubit decay rate so that we can understand their funda- 
mental physical origins as well as engineer improvements. 
Besides j K evaluated above, there are two additional con- 
tributions to the total damping rate 7 = j K + 7j_ + 7nr- 
Here 7^ is the decay rate into photon modes other than 
the cavity mode, and 7nr is the rate of other (possi- 
bly non-radiative) decays. Optical cavities are relatively 
open and 7_l is significant, but for ID microwave cav- 
ities, 7j_ is expected to be negligible (despite the very 
large transition dipole). For Rydberg atoms the two 
qubit states are both highly excited levels and 7nr rep- 
resents (radiative) decay out of the two-level subspace. 
For Cooper pair boxes, 7nr is completely unknown at the 
present time, but could have contributions from phonons, 
two-level systems in insulating | 2(| barriers and sub- 
strates, or thermally excited quasiparticles. 



For Cooper box qubits not inside a cavity, recent ex- 
periments |l8j have determined a relaxation time 1/7 = 
Ti ~ 1.3 fis despite the back action of continuous mea- 
surement by a SET electrometer. Vion et al. 01 found 
Ti ~ 1.84 [is (without measurement back action) for their 
charge-phase qubit. Thus in these experiments, if there 
are non-radiative decay channels, they are at most com- 
parable to the vacuum radiative decay rate (and may well 
be much less) estimated using Eq. (|18|) . Experiments 
with a cavity will present the qubit with a simple and 
well controlled electromagnetic environment, in which 
the radiative lifetime can be enhanced with detuning to 
l/7 re > 64 /xs, allowing 7nr to dominate and yielding 
valuable information about any non-radiative processes. 



VI. DISPERSIVE QND READOUT OF QUBIT 

In addition to lifetime enhancement, the dispersive 
regime is advantageous for read-out of the qubit. This 
can be realized by microwave irradiation of the cavity and 
then probing the transmitted or reflected photons |26J . 



A. Measurement Protocol 

A drive of frequency lo^ w on the resonator can be mod- 
eled by [HI 



0. 



(20) 



where e(t) is a measure of the drive the amplitude. In 
the dispersive limit, one expects from Fig. ^) peaks in 
the transmission spectrum at u> T — g 2 / A and fl + 2g 2 /A 
if the qubit is initially in its ground state. In a frame 
rotating at the drive frequency, the matrix elements for 
these transitions are respectively 



(T,0|/V|- n) 
(T.Olfl^Fhn) 



£ 

eg 

A' 



(21) 



In the large detuning case, the peak at £l + 2g 2 j A, corre- 
sponding approximatively to a qubit flip, is highly sup- 
pressed. 

The matrix element corresponding to a qubit flip from 
the excited state is also suppressed and, as shown in 
Fig. [5] depending on the qubit being in its ground or 
excited states, the transmission spectrum will present a 
peak of width k at lj t — g 2 /A or u> T + g 2 /A. With the 
parameters of Table this dispersive pull of the cavity 
frequency is ±g 2 / 'kA = ±2.5 line widths for a 10% de- 
tuning. Exact diagonalization Q shows that the pull is 
power dependent and decreases in magnitude for cavity 
photon numbers on the scale n — n CI i t = A 2 /Ag 2 . In the 
regime of non-linear response, single-atom optical bista- 
bility 14] can be expected when the drive frequency is off 
resonance at low power but on resonance at high power 

M. 
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FIG. 5: (color online). Transmission spectrum of the cavity, 
which is "pulled" by an amount ±g 2 /A = 2.5 x 10 -4 x ui t , 
depending on the state of the qubit (red for the excited state, 
blue for the ground state). To perform a measurement of the 
qubit, a pulse of microwave photons, at a probe frequency 
cj,j W = LU r or to T ±g 2 /A is sent through the cavity. Additional 
peaks near Q corresponding to qubit flips are suppressed by 



The state-dependent pull of the cavity frequency by 
the qubit can be used to entangle the state of the qubit 
with that of the photons transmitted or reflected by the 
resonator. For g 2 / kA > 1, as in Fig. EI the pull is greater 
than the line width and irradiating the cavity at one of 
the pulled frequencies lo t ± g /A, the transmission of the 
cavity will be close to unity for one state of the qubit and 
close to zero for the other |3(j ■ 

Choosing the drive to be instead at the bare cav- 
ity frequency uj r , the state of the qubit is encoded in 
the phase of the reflected and transmitted microwaves. 
An initial qubit state \x) = a|f) + evolves un- 

der microwave irradiation into the entangled state \ip) = 
a\1,0) +/3||,-0), where tan6> = 2g 2 /nA, and |±0) are 
(interaction representation) coherent states with the ap- 
propriate mean photon number and opposite phases. In 
the situation where g 2 / kA <C 1, this is the most appro- 
priate strategy. 

It is interesting to note that such an entangled state 
can be used to couple qubits in distant resonators and 
allow quantum communication |31| . Moreover, if an in- 
dependent measurement of the qubit state can be made, 
such states can be turned into photon Schrodinger cats 

m. 

To characterize these two measurement schemes cor- 
responding to two different choices of the drive fre- 
quency, we compute the average photon number inside 
the resonator n and the homodyne voltage on the 5017 
impedance at the output of the resonator. Since the 
power coupled to the outside of the resonator is P = 
(n)huj Y K/2 = (V ou t) 2 /R, the homodyne voltage can be 
expressed as (K>ut) = \ / Rhu r K(a + a))/2 and is propor- 
tional to the real part of the field inside the cavity. 

In the absence of dissipation, the time dependence of 
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FIG. 6: (color online). Results of numerical simulations us- 
ing the quantum state diffusion method. A microwave pulse 
of duration ~ 15/k and centered at the pulled frequency 
Lu r + g 2 /A drives the cavity, a) The occupation probability of 
the excited state (right axis), for the case in which the qubit 
is initially in the ground (blue) or excited (red) state and in- 
tracavity photon number (left axis), are shown as a function 
of time. Though the qubit states are temporarily coherently 
mixed during the pulse, the probability of real transitions is 
seen to be small. Depending on the qubit's state, the pulse is 
either on or away from the combined cavity-qubit resonance, 
and therefore is mostly transmitted or mostly reflected, b) 
The real component of the cavity electric field amplitude (left 
axis), and the transmitted voltage phasor (right axis) in the 
output transmission line, for the two possible initial qubit 
states. The parameters used for the simulation are presented 
in Tabled 



the field inside the cavity can be obtained in the Heisen- 
berg picture from Eqs. I|12|) and Q2U[I. This leads to a 
closed set of differential equations for a, a z and aa z which 
is easily solved. In the presence of dissipation however 
(i.e. performing the transformation 1)11(1 on H K and _ff 7 , 
and adding the resulting terms to Eqs. H12fl and l|2U|) ). 
the set is no longer closed and we resort to numerical 
stochastic wave function calculations p32| . 

Figures IHI and [7| show the numerical results for the two 
choices of drive frequency and using the parameters of 
Table |I] For these calculations, a pulse of duration ~ 
15/k with a hyperbolic tangent rise and fall, is used to 
excite the cavity. Fig. corresponds to a drive at the 
pulled frequency w r + g 2 /A. In Fig. E^,) the probability 
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FIG. 7: (color online). Same as Fig. [(J for the drive at the 
bare cavity frequency u) T . Depending on the qubit's state, the 
pulse is either above or below the combined cavity-qubit res- 
onance, and so is partly transmitted and reflected but with a 
large relative phase shift that can be detected with homodyne 
detection. In b), the opposing phase shifts cause a change in 
sign of the output, which can be measured with high signal- 
to-noise to realize a single-shot, QND measurement of the 
qubit. 



Pj to find the qubit in its excited state (right axis) is 
plotted as a function of time for the qubit initially in the 
ground (blue) or excited state (red). The dashed lines 
represent the corresponding number of photons in the 
cavity (left axis). Fig- El 3 ) shows, in a frame rotating at 
the drive frequency, the real part of the cavity electric 
field amplitude (left axis) and transmitted voltage phase 
(right axis) in the output transmission line, again for the 
two possible initial qubit states. These quantities are 
shown in Fig. for a drive at the bare frequency ui T . 

As expected, for the first choice of drive frequency, 
the information about the state of the qubit is mostly 
stored in the number of transmitted photons. When the 
drive is at the bare frequency however, there is very lit- 
tle information in the photon number, with most of the 
information being stored in the phase of the transmitted 
and reflected signal. This phase shift can be measured 
using standard heterodyne techniques. As also discussed 
in appendix [5] both approaches can serve as a high effi- 
ciency quantum non-demolition dispersive readout of the 
state of the qubit. 



B. Measurement Time and Backaction 

As seen from Eq. I|12(l . the back action of the disper- 
sive cQED measurement is due to quantum fluctuations 
of the number of photons n within the cavity. These fluc- 
tuations cause variations in the ac Stark shift (g 2 / A)na z 
that in turn dephase the qubit. It is useful to compute 
the corresponding dephasing rate and compare it with 
the measurement rate, i.e. the rate at which information 
about the state of the qubit can be acquired. 

To determine the dephasing rate, we assume that the 
cavity is driven at the bare cavity resonance frequency 
and that the pull of the resonance is small compared to 
the line width k. The relative phase accumulated be- 
tween the ground and excited states of the qubit is 

^(t) = 2^- f dt'n{t') (22) 
A Jo 

which yield a mean phase advance (ip) = 29qN with 
#o = 2g 2 /kA and N = Kfvt/2 the total number of trans- 
mitted photons |14| . For weak coupling, the dephasing 
time will greatly exceed 1/k and, in the long time limit, 
the noise in ip induced by the ac Stark shift will be gaus- 
sian. Dephasing can then be evaluated by computing the 
long time decay of the correlator 

~ e -K 2 ^) 2 IS /o «Mi«Ma<»(tiMta)> (23) 

To evaluate this correlator in the presence of a 
continuous- wave (CW) drive on the cavity, we first per- 
form a canonical transformation on the cavity operators 
<j(t) by writing them in terms of a classical and a 
quantum part 

a(t) =a(t)+d(t). (24) 

Under this transformation, the coherent state obeying 
a | a) = a \ a), is simply the vacuum for the operator d. It 
is then easy to verify that 

((n(t) - ft)(n(0) - ft)) = a 2 (d(t)d^(0)) = ne"* 1 ' 1 . (25) 

It is interesting to note that the factor of 1/2 in the ex- 
ponent is due to the presence of the coherent drive. If 
the resonator is not driven, the photon number correlator 
rather decays at a rate k. Using this result in ill' Ml yields 
the dephasing rate 

r„=40g£n. (26) 

Since the rate of transmission on resonance is ren/2, this 
means that the dephasing per transmitted photon is W 2 . 

To compare this result to the measurement time Tineas; 
we imagine a homodyne measurement to determine the 
transmitted phase. Standard analysis of such an interfer- 
ometric set up |l4j shows that the minimum phase change 
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which can be resolved using N photons is 59 = l/\/~N. 
Hence the measurement time to resolve the phase change 
59 = 29 is 



which yields 

T m rv = 1- (28) 

This exceeds the quantum limit |33j T m T v = 1/2 by a fac- 
tor of 2. Equivalently, in the language of Ref. [34| (which 
uses a definition of the measurement time twice as large 
as that above) the efficiency ratio is x = 1/PmIV) = 0.5. 

The failure to reach the quantum limit can be traced 
[35| to the fact that that the coupling of the photons 
to the qubit is not adiabatic. A small fraction R m 9q 
of the photons incident on the resonator are reflected 
rather than transmitted. Because the phase shift of the 
reflected wave differs by ir between the two states of 
the qubit, it turns out that, despite its weak intensity, 
the reflected wave contains precisely the same amount of 
information about the state of the qubit as the transmit- 
ted wave which is more intense but has a smaller phase 
shift. In the language of Ref. [34|. this 'wasted' infor- 
mation accounts for the excess dephasing relative to the 
measurement rate. By measuring also the phase shift of 
the reflected photons, it could be possible to reach the 
quantum limit. 

Another form of possible back action is mixing tran- 
sitions between the two qubit states induced by the mi- 
crowaves. First, as seen from Fig. Et) and [7^), increasing 
the average number of photons in the cavity induces mix- 
ing. This is simply caused by dressing of the qubit by the 
cavity photons. Using the dressed states @ and ©, the 
level of this coherent mixing can be estimated as 

P U = i(±^|l±^|±^) (29) 

= \{l±- = ) (30) 

Exciting the cavity to n = n cr ;t, yields Pj_ ~ 0.85. As 
is clear from the numerical results, this process is com- 
pletely reversible and does not lead to errors in the read- 
out. 

The drive can also lead to real transitions between the 
qubit states. However, since the coupling is so strong, 
large detuning A = 0.1 cj r can be chosen, making the 
mixing rate limited not by the frequency spread of the 
drive pulse, but rather by the width of the qubit ex- 
cited state itself. The rate of driving the qubit from 
ground to excited state when n photons are in the cavity 
is R ~ n(g/ A) 2r y. If the measurement pulse excites the 
cavity to n — n cr it, we see that the excitation rate is still 
only 1/4 of the relaxation rate. As a result, the main 
limitation on the fidelity of this QND readout is the de- 
cay of the excited state of the qubit during the course of 



the readout. This occurs (for small 7) with probability 
-frciax ~ 7imeas ~ 15 x j / k ~ 3.75 % and leads to a small 
error P crr ~ 57/K ~ 1.5% in the measurement, where 
we have taken 7 = j K . As confirmed by the numerical 
calculations of Fig. [5] and this dispersive measurement 
is therefore highly non-demolition. 



C. Signal-to-Noise 

For homodyne detection in the case where the cavity 
pull g 2 /A«; is larger than one, the signal-to-noise ratio 
(SNR) is given by the ratio of the number of photons 
n s ig = nnAt/2 accumulated over an integration period 
At, divided by the detector noise n amp = fcsT/v/?iu; r . 
Assuming the integration time to be limited by the 
qubit 's decay time 1/7 and exciting the cavity to a max- 
imal amplitude n cr ;t = 100 ~ ?T- amp , we obtain SNR = 
("crit/^amp)(K/27). If the qubit lifetime is longer than 
a few cavity decay times (1/k = 160ns), this SNR can 
be very large. In the most optimistic situation where 
7 = 7 K , the signal-to- noise ratio is SNR=200. 

When taking into account the fact that the qubit has 
a finite probability to decay during the measurement, a 
better strategy than integrating the signal for a long time 
is to take advantage of the large SNR to measure quickly. 
Simulations have shown that in the situation where 7 = 
7 K , the optimum integration time is roughly 15 cavity 
lifetimes. This is the pulse length used for the stochastic 
numerical simulations shown above. The readout fidelity, 
including the effects of this stochastic decay, and related 
figures of merit of the single-shot high efficiency QND 
readout are summarized in Table ITT1 

This scheme has other interesting features that are 
worth mentioning here. First, since nearly all the energy 
used in this dispersive measurement scheme is dissipated 
in the remote terminations of the input and output trans- 
mission lines, it has the practical advantage of avoiding 
quasiparticle generation in the qubit. 

Another key feature of the cavity QED readout is that 
it lends itself naturally to operation of the box at the 
charge degeneracy point (N g = 1/2), where it has been 
shown that T2 can be enormously enhanced 01 because 
the energy splitting has an extremum with respect to 
gate voltage and isolation of the qubit from 1/f dephas- 
ing is optimal. The derivative of the energy splitting 
with respect to gate voltage is the charge difference in the 
two qubit states. At the degeneracy point this derivative 
vanishes and the environment cannot distinguish the two 
states and thus cannot dephase the qubit. This also im- 
plies that a charge measurement cannot be used to deter- 
mine the state of the system 0,0- While the first deriva- 
tive of the energy splitting with respect to gate voltage 
vanishes at the degeneracy point, the second derivative, 
corresponding to the difference in charge polarizability 
of the two quantum states, is maximal. One can think 
of the qubit as a non-linear quantum system having a 
state-dependent capacitance (or in general, an admit- 
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parameter 


symbol 


ID circuit 


dimensionless cavity pull 


g 2 / K A 


2.5 


cavity-enhanced lifetime 


j- 1 = (A/gfn- 1 


64 ^s 


readout SNR 


SNR = (7i cr it/n amp )«:/27 


200 (6) 


readout error 


P crr ~5x 7/k 


1.5 % (14%) 


1 bit operation time 


T„ > 1/A 


> 0.16 ns 


entanglement time 


*x/?Swap - 7rA/4g 2 


~ 0.05 (js 


2 bit operations 


iVop - l/[7ViSWAp] 


> 1200 (40) 



TABLE II: Figures of merit for readout and multi-qubit entanglement of superconducting qubits using dispersive (off-resonant) 
coupling to a ID transmission line resonator. The same parameters as Table 1, and a detuning of the Cooper pair box from 
the resonator of 10% (A = 0.1 u> T ), are assumed. Quantities involving the qubit decay 7 are computed both for the theoretical 
lower bound 7 = 7 K for spontaneous emission via the cavity, and (in parentheses) for the current experimental upper bound 
1/7 > 2(js. Though the signal-to-noise of the readout is very high in either case, the estimate of the readout error rate 
is dominated by the probability of qubit relaxation during the measurement, which has a duration of a few cavity lifetimes 
(~ 1 - 10 k" 1 ). If the qubit non-radiative decay is low, both high efficiency readout and more than 10 3 two-bit operations 
could be attained. 



tance) which changes sign between the ground and ex- 
cited states [3(| . It is this change in polarizability which 
is measured in the dispersive QND measurement. 

In contrast, standard charge measurement schemes 
[lsL I]}?! require moving away from the optimal point. 
Simmonds et al. |2(j have recently raised the possibil- 
ity that there are numerous parasitic environmental res- 
onances which can relax the qubit when its frequency 
f2 is changed during the course of moving the operat- 
ing point. The dispersive cQED measurement is there- 
fore highly advantageous since it operates best at the 
charge degeneracy point. In general, such a measurement 
of an ac property of the qubit is strongly desirable in 
the usual case where dephasing is dominated by low fre- 
quency (1/f) noise. Notice also that the proposed quan- 
tum non-demolition measurement would be the inverse 
of the atomic microwave cQED measurement in which 
the state of the photon field is inferred non-destructively 
from the phase shift in the state of atoms sent through 
the cavity 0]. 



VII. COHERENT CONTROL 

While microwave irradiation of the cavity at its reso- 
nance frequency constitutes a measurement, irradiation 
close to the qubit 's frequency can be used to coherently 
control the state of the qubit. In the former case, the 
phase shift of the transmitted wave is strongly dependent 
on the state of the qubit and hence the photons become 
entangled with the qubit, as shown in Fig.|S| In the latter 
case however, driving is not a measurement because, for 
large detuning, the photons are largely reflected with a 
phase shift which is independent of the state of the qubit. 
There is therefore little entanglement between the field 
and the qubit in this situation and the rotation fidelity 
is high. 

To model the effect of the drive on the qubit, we 



add the microwave drive of Eq. I|2U|I to the Jaynes- 
Cumming Hamiltonian Q and apply the transformation 
(|llfl (again neglecting damping) to obtain the following 
effective one-qubit Hamiltonian 
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Lo^ w )a)a + He(t)((r + a), 



(31) 



in a frame rotating at the drive frequency lo^ w . Choos- 
ing uj^ w = Q + (2n + l)g 2 /A, Hi q generates rotations of 
the qubit about the x axis with Rabi frequency gs/A. 
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FIG. 8: (color online). Phase shift of the cavity field for the 
two states of the qubit as a function of detuning between the 
driving and resonator frequencies. Obtained from the steady- 
state solution of the equation of motion for a(t) while only 
taking into account damping on the cavity and using the pa- 
rameters of Table H] Read-out of the qubit is realized at, or 
close to, zero detuning between the drive and resonator fre- 
quencies where the dependence of the phase shift on the qubit 
state is largest. Coherent manipulations of the qubit are real- 
ized close to the qubit frequency which is 10% detuned from 
the cavity (not shown on this scale). At such large detunings, 
these is little dependence of the phase shift on the qubit 's 
state. 
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Different drive frequencies can be chosen to realize ro- 
tations around arbitrary axes in the x-z plane. In par- 
ticular, choosing lj^ w = Q + [2n + \)g 2 / A — 2ge/A and 
t = 7rA/2v / 2ge generates the Hadamard transformation 
H. Since 7ia x H = a z , these two choices of frequency are 
sufficient to realize any 1-qubit logical operation. 

Assuming that we can take full advantage of lifetime 
enhancement inside the cavity (i.e. that 7 = j K ), the 
number of it rotations about the x axis which can be car- 
ried out is — 2eA/irgK ~ 10 5 e for the experimental 
parameters assumed in Table Q] For large e, the choice 
of drive frequency must take into account the power de- 
pendence of the cavity frequency pulling. 

Numerical simulation shown in Fig. [5] confirms this 
simple picture and that single-bit rotations can be per- 
formed with very high fidelity. It is interesting to note 
that since detuning between the resonator and the drive 
is large, the cavity is only virtually populated, with an 
average photon number n « e 2 / A 2 ~ 0.1. Virtual pop- 
ulation and depopulation of the cavity can be realized 
much faster than the cavity lifetime 1/k and, as a result, 
the qubit feels the effect of the drive rapidly after the 
drive has been turned on. The limit on the speed of turn 
on and off of the drive is set by the detuning A. If the 
drive is turned on faster than 1/A, the frequency spread 
of the drive is such that part of the drive's photons will 
pick up phase information (see Fig. |SJ and dephase the 
qubit. As a result, for large detuning, this approach leads 
to a fast and accurate way to coherently control the state 
of the qubit. 

To model the effect of the drive on the resonator an 
alternative model is to use the cavity-modified Maxwell- 
Bloch equations [2^]. As expected, numerical integra- 
tion of the Maxwell-Bloch equations reproduce very well 
the stochastic numerical results when the drive is at the 
qubit 's frequency but do not reproduce these numerical 
results when the drive is close to the bare resonator fre- 
quency (Fig. and[7J, i.e. when entanglement between 
the qubit and the photons cannot be neglected. 



VIII. RESONATOR AS QUANTUM BUS: 
ENTANGLEMENT OF MULTIPLE QUBITS 

The transmission-line resonator has the advantage that 
it should be possible to place multiple qubits along its 
length (~ 1 cm) and entangle them together, which is an 
essential requirement for quantum computation. For the 
case of two qubits, they can be placed closer to the ends of 
the resonator but still well isolated from the environment 
and can be separately dc biased by capacitive coupling to 
the left and right center conductors of the transmission 
line. Additional qubits would have to have separate gate 
bias lines installed. 

For the pair of qubits labeled i and j, both coupled 
with strength g to the cavity and detuned from the res- 
onator but in resonance with each other, the transfor- 
mation (|ll|l yields the effective two-qubit Hamiltonian 
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FIG. 9: (color online). Numerical stochastic wave function 
simulation showing coherent control of a qubit by microwave 
irradiation of the cavity at the ac-Stark and Lamb shifted 
qubit frequency. The qubit is first left to evolve freely for 
about 40ns. The drive is turned on for t = TiiA/2ge ~ 115ns, 
corresponding to 7n pulses, and then turned off. Since the 
drive is tuned far away from the cavity, the cavity photon 
number is small even for the moderately large drive amplitude 
e = 0.03oj r used here. 
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In addition to ac-Stark and Lamb shifts, the last term 
couples the qubits thought virtual excitations of the res- 
onator. 

In a frame rotating at the qubit's frequency Q, Hi q 
generates the evolution 
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where l r , is the identity operator in the resonator space. 
Up to phase factors, this corresponds at t = irA/ig 2 ~ 
50 ns to a %/ iSWAP logical operation. Up to one-qubit 
gates, this operation is equivalent to the controlled-NOT. 
Together with one-qubit gates, the interaction Hiq is 
therefore sufficient for universal quantum computation 
[4fT |. Assuming again that we can take full advantage 
of the lifetime enhancement inside the cavity, the num- 
ber of VzSWAP operations which can be carried out is 
= 4A/7TK ~ 1200 for the parameters assumed above. 
This can be further improved if the qubit's non-radiative 
decay is sufficiently small, and higher Q cavities are em- 
ployed. 
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When the qubits are detuned from each other, the off- 
diagonal coupling provided by H 2q is only weakly effec- 
tive and the coupling is for all practical purposes turned 
off. Two-qubit logical gates in this setup can therefore be 
controlled by individually tuning the qubits. Moreover, 
single-qubit and two-qubit logical operations on different 
qubits and pairs of qubits can both be realized simultane- 
ously, a requirement to reach presently known thresholds 
for fault-tolerant quantum computation |4l| . 

It is interesting to point out that the dispersive QND 
readout presented in section lVll mav be able to determine 
the state of multiple qubits in a single shot without the 
need for additional signal ports. For example, for the 
case of two qubits with different detunings, the cavity 
pull will take four different values ±g , f/Ai±g , |/A2 allow- 
ing single-shot readout of the coupled system. This can 
in principle be extended to N qubits provided that the 
range of individual cavity pulls can be made large enough 
to distinguish all the combinations. Alternatively, one 
could read them out in small groups at the expense of 
having to electrically vary the detuning of each group to 
bring them into strong coupling with the resonator. 



IX. ENCODED UNIVERSALITY AND 
DECOHERENCE-FREE SUBSPACE 

Universal quantum computation can also be realized 
in this architecture under the encoding C = {|ti) , lit)} 
by controlling only the qubit's detuning and, therefore, 
by turning on and off the interaction term in H^q |42|. 

An alternative encoded two-qubit logical operation 
to the one suggested in Ref. |42| can be realized here 
by tuning the four qubits forming the pair of encoded 
qubits in resonance for a time t — irA/3g 2 . The re- 
sulting effective evolution operator can be written as 
U 2q = cxp [-i(irA/3g 2 )a 

xi&x2 j where a X i is a Pauli op- 
erator acting on the i th encoded qubit. Together with 
encoded one-qubit operations, U?, q is sufficient for uni- 
versal quantum computation using the encoding C. 

We point out that the subspace £ is a decoherence-free 
subspace with respect to global dephasing and use of 
this encoding will provide some protection against noise. 
The application of L^g on the encoded subspace C how- 
ever causes temporary leakage out of this protected sub- 
space. This is also the case with the approach of Ref. ^2 • 
In the present situation however, since the Hamiltonian 
generating L^q commutes with the generator of global 
dephasing, this temporary excursion out of the protected 
subspace does not induce noise on the encoded qubit. 

X. SUMMARY AND CONCLUSIONS 

In summary, we propose that the combination of 
one-dimensional superconducting transmission line res- 
onators, which confine their zero point energy to ex- 
tremely small volumes, and superconducting charge 



qubits, which are electrically controllable qubits with 
large electric dipole moments, constitutes an interest- 
ing system to access the strong-coupling regime of cavity 
quantum electrodynamics. This combined system is an 
advantageous architecture for the coherent control, en- 
tanglement, and readout of quantum bits for quantum 
computation and communication. Among the practical 
benefits of this approach are the ability to suppress radia- 
tive decay of the qubit while still allowing one-bit opera- 
tions, a simple and minimally disruptive method for read- 
out of single and multiple qubits, and the ability to gener- 
ate tunable two-qubit entanglement over centimeter-scale 
distances. We also note that in the structures described 
here, the emission or absorption of a single photon by the 
qubit is tagged by a sudden large change in the resonator 
transmission properties 29] making them potentially use- 
ful as single photon sources and detectors. 
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APPENDIX A: QUANTIZATION OF THE ID 
TRANSMISSION LINE RESONATOR 

A transmission line of length L, whose cross section 
dimension is much less then the wavelength of the trans- 
mitted signal can be approximated by a 1-D model. For 
relatively low frequencies it is well described by an infinite 
series of inductors with each node capacitively connected 
to ground, as shown in Fig. [21 Denoting the inductance 
per unit length I and the capacitance per unit length c, 
the Lagrangian of the circuit is 

where j(x, t) and q(x, t) are the local current and charge 
density, respectively. We have ignored for the moment 
the two semi-infinite transmission lines capacitively cou- 
pled to the resonator. Defining the variable 0(x,i) 

6{x,t) = [ dx'q(x',t), (A2) 

J-L/2 
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the Lagrangian can be rewritten as 



C 



L/2 / 1 ■ 1 

dx ( V_ J_(V0) 2 

-L/2 V 2 



2c 



(A3) 



The corresponding Euler-Lagrange equation is a wave 
equation with the speed v = Wl/lc. Using the boundary 
conditions due to charge neutrality 



6(-L/2,t) = 0(L/2,t) = 0, 



(A4) 



we obtain 



0(x,t) 



_ feo, cutoff , 

i V—v , . K TTX 

4>k a (t) cos 



^e, cutoff 



7 X! (*) sin 



fc c =2 



L 

k c TTX 

TT' 



(A5) 



for odd and even modes, respectively. For finite length 
L, the transmission line acts as a resonator with resonant 
frequencies lju — knv/L. The cutoff is determined by the 
fact that the resonator is not strictly one dimensional. 

Using the normal mode expansion (|A5|) in (|A3|) . one 
obtains, after spatial integration, the Lagrangian in the 
form of a set of harmonic oscillators 



£ = 2^2** ~ Yc It 



(A6) 



Promoting the variable 0fc and its canonically conju- 
gated momentum ix^ = l(f>). to conjugate operators and 
introducing the boson creation and annihilation opera- 
tors at and a./, satisfying [cifc,a£,] = S^i, we obtain the 
usual relations diagonalizing the Hamiltonian obtained 
from the Lagrangian (|A6|) 



.(*) = 



fthJhC L 



(a k (t)+al(t)) (A7) 



2 ten 

**(*) = -iY^(°*(*)-4(*)). (A8) 

From these relations, the voltage on the resonator can be 
expressed as 



V(x,t) 



1 d6{x,t) 
c dx 



£ 

fco = l 



TT 



k TTX 



(A9) 

K(*) + 4,(*)] 



fcc = l 



In the presence of the two semi-infinite transmission 
lines coupled to the resonator, the Lagrangian 1A3|I and 
the boundary conditions (|A4|) are modified to take into 
account the voltage drop on the coupling capacitors Cq. 
Assuming no spatial extent for the capacitors Co, the 
problem is still solvable analytically. Due to this cou- 
pling, the wavefunction can now extend outside of the 
central segment which causes a slight red-shift, of order 
Cq/Lc, of the cavity resonant frequency. 

As shown in Fig. [21 we assume the qubit to be fabri- 
cated at the center of the resonator. As a result, at low 
temperatures, the qubit is coupled to the mode k = 2 
of the resonator, which as an anti-node of the voltage 
in its center. The rms voltage between the cen ter con- 
ductor and the ground plane is then V^ s = \J huj v /cL 
with uj i = uj 2 and the voltage felt by the qubit is 

i/(o, t) = Kms( a 2(*) + <4(*))- In tne main bod y of tnis 

paper, we work only with this second harmonic and drop 
the mode index on the resonator operators. 



APPENDIX B: QUANTUM NON-DEMOLITION 
MEASUREMENTS 



Read-out of a qubit can lead to both mixing and de- 
phasing [23l l33l | . While dephasing is unavoidable, mixing 
of the measured observable can be eliminated in a QND 
measurement by choosing the qubit-measurement appa- 
ratus interaction such that the measured observable is a 
constant of motion. In that situation, the measurement- 
induced mixing is rather introduced in the operator con- 
jugate to the operator being measured. 

In the situation of interest in this paper, the op- 
erator being probed is a z and, from Eq. i|12[l . the 
qubit-measurement apparatus interaction Hamiltonian is 
-Hint = (5 2 /A)cr z a 1 'a, such that [a z ,Hi n t] = 0. For a z to 
be a constant of motion also requires that it commutes 
with the qubit Hamiltonian. This condition is also satis- 
fied in Eq. ifTl . 

That the measured observable is a constant of motion 
implies that repeated observations will yield the same 
result. This allows for the measurement result to reach 
arbitrary large accuracy by accumulating signal. In prac- 
tice however, there are always environmental dissipation 
mechanisms acting on the qubit independently of the 
read-out. Even in a QND situation, these will lead to 
a finite mixing rate 1/T\ of the qubit in the course of the 
measurement. Hence, high fidelity can only be achieved 
by a strong measurement completed in a time T m <C T±. 
This simple point is not as widely appreciated as it should 
be. 
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